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Abstract 

This is the first of a series of two papers on the perturbative computation of 
the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Our final aim is 
to eventually compare with results from lattice simulations in order to enlight the 
genuinely non-perturbative content of the latter. By means of Numerical Stochastic 
Perturbation Theory we compute the ghost propagator in Landau gauge up to three 
loops. We present results in the infinite volume and a — > limits, based on a general 
strategy that we discuss in detail. 



1 Introduction 

The identification of confinement mechanisms is a challenge for Lattice QCD [H El IB] ■ This 
requires to single out observables which both capture the relevant physical content and 
are in the end suitable for numerical simulations, so that a non-perturbative investigation 
becomes possible. The gluon and the ghost propagators with their momentum dependence 
are believed to encode specific features and mechanisms of a confining theory [2] , which can 
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be studied once a gauge is fixed. Apart from the non-unique asymptotic infrared behavior, 
which - if calculated in lattice simulations - realizes the so-called decoupling solution 
instead of the scaling solution (see Refs. [U EJ [6]) of the Schwinger-Dyson and Functional 
Renormalization Group [7J approaches 0, in these functions the onset of non-perturbative 
effects is generally agreed to be associated with the following two aspects: the coupling to 
local condensates [13] and to localized vacuum excitations like vortices [HJ [TJ] . In order to 
effectively face the genuine non-perturbative effects, it is desirable to have for reference an 
understanding of the perturbative behavior of these gauge-dependent functions in higher 
than one-loop order. 

Such a program cannot evade a well-known practical problem: Lattice Perturbation 
Theory (LPT) is much more involved compared to its continuum QCD counterpart. The 
complexity of the diagrammatic approaches rapidly increases beyond the one-loop approx- 
imation [16]. By now, only a limited number of results up to two-loop accuracy has been 
obtained. As an effective alternative to diagrammatic LPT a numerical tool has been pro- 
posed by the Parma group in [T7] . The method is a quite direct numerical implementation 
of Stochastic Perturbation Theory, which in turn comes as one of the two main applications 
of the ideas of Stochastic Quantization, the other being the use of Langevin dynamics as 
an alternative to standard Monte Carlo techniques [TBI [19] . This numerical tool for LPT 
goes under the name of Numerical Stochastic Perturbation Theory (NSPT) and has been 
reviewed in |20|. 



In a joined effort we started a new field of application for NSPT based on our respective 
interests in LPT (Leipzig), NSPT (Parma) and Monte Carlo studies of gluon and ghost 
propagators (Berlin). First results have been presented at Lattice 2007 [2U [22], Lattice 
2008 [23J, Confinements [21] and at Lattice 2009 [25]. The present paper is the first of a 
series of two. While this paper is concentrated on the ghost propagator, the forthcoming 
one will present the results for the gluon propagator. 

One of our main goals is to provide results for the comparison of higher-loop calculations 
of perturbative contributions with the results of lattice Monte Carlo calculations of the 
main propagators in QCD on finite size lattices. The other goal consists in determining 
higher loop contributions to the logarithmically divergent perturbative ghost propagator in 
momentum space including non-logarithmic contributions to good precision. The present 
work also aims at a further check on the precision of the method by comparison to known 
one-loop results from LPT. In order to proceed to two- and higher-loop results in the 
infinite volume limit, peculiar emphasis will be put on the control of finite size effects. 
We want to stress that our results have been obtained from two completely independent 
sources: while strategies have been extensively discussed, the actual implementation of 
Parma and Leipzig NSPT codes are completely independent of each other, thus providing 
a stronger check on numerical results. 

The paper is organized as follows. In Sect. [2] we present selected elements of NSPT 

1 The latter solution is preferred in view of the Gribov-Zwanziger j8j [9] and Kugo-Ojima [10l QJ] 
confinement scenarios. The conflict with lattice calculations has given rise to a controversy about the 
(in)dispcnsability of having a BRST-invariant lattice gauge-fixing prescription. It is argued to be possi- 
ble [12] although still lacking. This interesting development is probably irrelevant for the aspects discussed 
below. 
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which are relevant for the case at hand, in particular the issue of performing a NSPT 
computation in a fixed gauge, i.e. the minimal Landau gauge. Sect. [3] focuses on the 
observables we will be concerned with, while Sect. H] summarizes known analytic results 
for the perturbative ghost propagator up to three-loop order and their relation to the 
unknown non-logarithmic coefficients we want to determine. In the following section we 
present the technical details of practicing NSPT. Finally we collect our results in Sect. [H] 
and briefly discuss a comparison to Monte Carlo data. 



2 Landau gauge fixing in NSPT 

The basis of Stochastic Quantization is the (stochastic) time evolution of fields driven by 
the Langevin equation. In the case of a SU(N) lattice gauge theory the Langevin equation 
for the SU(N) matrix U Xjfl located at the link between lattice sites x and x + jl reads 

^U x ^{t-r 1 ) = i{V x ^S G [U}-r lx ^t)) U x ^(t; V ) , (1) 

t is the (stochastic) Langevin time. The so-called drift term is given by the equations 
of motion derived from the Euclidean gauge action Sg[U]: they are written in terms of 
the differential operator V XiAt = T a V% ^ which is a left Lie derivative. The T a (a = 
1, . . . , N 2 — 1) are the generators in the fundamental representation. In our case we will 
have N = 3 and the gauge action Sg[U] is the Wilson plaquette action. 

The key issue of Stochastic Quantization is that in the limit of large t expectation 
values with respect to the noise converge to expectation values with respect to the QFT 
functional integral, i.e. the distribution of the gauge fields converges to the Gibbs measure 
P[U] oc exp(— Sa[U]). In any numerical approach (and this applies to our case), however, 
the Langevin equation is not solved for continuous time. One needs to discretize time as 
a sequence t = ne, with running step number n. Therefore, in order to extract correct 
physical information, one needs to perform the double extrapolation n — > oo and e — > 0. 
For the numerical solution of the Langevin equation we adhere to a peculiar version of 
the Euler scheme that guarantees all the link matrices U Xjfl G SU (3) to stay in the group 
manifold: 

U Xtft [n + l;ri) = exp(i F x>(t [U, rj]) U x ^(n; rj) , (2) 
F x jU,r]} = eV Xt(t S G [U] + VZriw (3) 

In matrix form and at finite e, the components of rj are subject to the mutual correlations 
(i, k, I, m — 1, . . . , 3) 



(Vi,k( Z ) Vl,m( W ))ri 



S zw . (4) 



z and w (integer) denote here Langevin time steps. The noise is white, meaning that 
different time steps are uncorrelated. 

The usual path from Stochastic Quantization to Stochastic Perturbation Theory would 
go through trading the differential form of Langevin equation for its integral equation 
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counterpart in a Green function formalism, which is easy to solve in a perturbative scheme. 
NSPT remains instead in a differential formulation. For the case of Lattice Gauge Theories 
we rewrite each link matrix as an expansion in the bare coupling constant go- Since 
ft = 2N/gQ, the expansion reads 

U^fati^l + ^P^UWfar,). (5) 

Provided one rescales the time step to e — fte, the expansion converts the Langevin equation 
(J2J) into a system of simultaneous updates in terms of the expansion coefficients (jSJ) of U XtfJ , 
and of similar expansion coefficients for the force F Xjti in Eq. (j3J). While the random noise 77 
enters only the lowest order equation, higher orders are stochastic by the noise propagating 
from lower order to higher order terms. The system can be truncated according to the 
maximal order of the gauge field one is interested in. 

In a notation which is common in Lattice Perturbation Theory we can relate the (an- 
tihermitean) vector potential A x +p,/2 >fi to the link matrices U XIJt via 

= log U XtfM . (6) 

An expansion similar to (jSJ) exists for that potential taking values in the algebra su(3), 

; /^*2W*J v) ■ (7) 

i>i 

We observe that unitarity of U XilM has a counterpart in the antihermiticity and tracelessness 
of all orders of A x+ ^/2 :fl : 

/t(0t -_4(0 Tr 4® -n (R\ 

As we will see, the definition of the ghost propagator will in a sense be based on the 
expansion in the algebra. 

We finally remark that whenever we speak about contributions of some order to an 
observable this has to be understood as an expansion 

(o)^J2p~ 1/2 (° {1) ) > ( 9 ) 

and constructed by comparison to coefficients of ft~ l l 2 . 

We now come to the main point of how to perform a gauge-fixed NSPT computation, 
in the case at hand in Landau gauge. A key issue of the standard (in other words, non- 
perturbative) Langevin equation for LGT is the non-necessity of fixing a gauge. In short: 
the drift term in Langevin equation is equal to the lattice equations of motion (put equal 
zero there). As a consequence, longitudinal degrees of freedom undergo an unbounded 
diffusion, but the group manifold is compact and hence no divergence shows up. Such a 
situation is no longer true in NSPT. This is quite obvious in the algebra: the diffusion of 
the longitudinal component of the fields is unbounded and hence their norms diverge in 
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stochastic time. We stress that gauge-invariant quantities are in principle not affected by 
these divergences, but this is no advantage: since norms are diverging one runs eventually 
into trouble by lack of accuracy. We also remark that this is true regardless of the choice 
of either expansion (17)) or (151): perturbation theory implies decompactification in any case. 

A solution comes along the lines suggested in : fixing (to some precision) the gauge 
is going to fix the problem and this goes under the name of Stochastic Gauge Fixing. We 
adhere here to the recipe of [20J. The Langevin evolution steps are interwoven with gauge 
transformations 

U^ = G{x)U x ^G\x + fi). (10) 

The main requirement for G is quite obvious: gauge transformations should be expanded 
in the coupling as well and the leading order should not affect the background vacuum 
around which the solution to Langevin equation is expanded (the trivial vacuum, in the 
case at hand). A convenient solution comes with the choice 



G(x)=exp I -a ] (*) ■ ( n ) 




In general the action of the left/right lattice derivative d^ R on a function f(x) defined 
on site x is given by 

di /R f{x)=±-(J{x)-f{xTp)) • (12) 

The action of the left lattice derivative on the vector potential defined on midpoints of the 
links is defined as 



The resulting divergence is mapped to site x. We will come back to this choice: the 
iteration of this gauge transformation would drive the field configuration to the Landau 
gauge. 

With respect to Stochastic Gauge Fixing, it suffices here to point out that one step of 
f fTUj) using (jlip interwoven with the Langevin step (T2]) is enough to keep fluctuations under 
control provided that a is chosen of order a ~ e. The rationale for this is the essence of 
Stochastic Gauge Fixing: any new drift term in the Langevin equation which amounts to 
a local gauge transformation does not affect the asymptotic distribution of gauge invariant 
quantities, while it can provide a restoring force in the non-gauge-invariant sector. The 
differential (finite difference, actually) operator present in (fTTj) makes it clear that zero 
modes are completely unaffected up to this point. Their divergences are still present. We 
take the simplest prescription of subtracting zero modes at every order. 

One can actually do more than this and perform a gauge-fixed NSPT computation in 
the case of covariant gauges. The interested reader is referred to [27]. The key issue is 
that one can rewrite the functional integral in the Faddeev-Popov formalism, which merely 
results in a different form of the gauge action (an effective action, one would say). Without 
entering into details, in this context it is worth making two points. First, there is no need 
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for the introduction of ghost fields: one simply faces an extra contribution to the action 
in the form Sfp = TrlogM, M being the Faddeev- Popov operator we will be concerned 
with in the following. Second, for the Landau gauge the procedure is of little help, since 
the other extra piece in the action is the gauge fixing action S g f = l/(2£) ^2 x Tr(d ■ A) 2 , 
with £ being the covariant gauge parameter, which is zero in this case. 

The failure of the general covariant gauge NSPT procedure has anyway a counterpart: 
as a limit case £ — > of covariant gauges, the Landau gauge 



(14) 



is related to an extremum of the norm of the field, and there is a well-known iterative 
procedure to reach such an extremum. The basic step is fTTTl) . i.e. the gauge transformation 
that we also use for Stochastic Gauge Fixing. The iteration stops when, for all orders Z, 
the summed-up quadratic violations of ( jLfl) vanish within double precision (V is the lattice 
volume), 



1-1 



[x 




0. 



(15) 



Iterating (TTOT) with (fTTj) is actually not very efficient, since it suffers from critical slowing 
down due to the low modes, which becomes more and more severe as the lattice size 
increases. 

An effective solution to this problem is also well-known and goes through Fourier ac- 
celeration as introduced in [28J. In NSPT one makes use of the expanded variant of the 
gauge transformation 



G{x) = exp { F~ 



a—-F 



5>H 



x 



(p) 



(16) 



in order to achieve the minimal Landau gauge. F denotes the Fourier transformation 
from discrete x to momentum p space and F~ l its inverse. The gauge transformation ([16]) 
applied locally is actually not local in the sense of ([TT]) . The driving force, the divergence of 
the gauge potential, is non-locally smeared with a low-pass filter that attenuates the higher 
Fourier components of d^A® by a factor p 2 m& Jp 2 ■ Here p 2 is a non-zero eigenvalue of the 
free field lattice Laplacian, and its maximal value on a periodic lattice in all directions is 
a2 Pmax = 16. While we notice that the volume- independent optimal parameter a for the 
case of periodic boundary conditions in all directions is found for this perturbative Landau 
gauge fixing to be a = l/(a 2 p^ ax ), we also point out that Fourier acceleration also works if 
one takes the plain lattice momenta (which in the following we will reference as p M ) instead 
of ?v 

This is a good point to pin down some notations for lattice momenta. In the following 
we will reserve the latin letter k to denote k^, the integer 4-tuples (ki, k 2 , k 3 , /c 4 ), while 
plain lattice momenta will be denoted by p^{k^) 

PM = ^. (17) 
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Here L M is the lattice size in direction \i and the reference to lattice spacing a can be made 
explicit: = N^a. In our case lattices will be symmetric, i.e. = L and = N. 
According to the two most popular conventions, G (— or k^ G [0, — 1]. 
The eigenvalues of the free Laplacian are 



where for each direction separately 



p„( W = W^=?sin(^). (19) 



a J a V 2 

3 The lattice ghost propagator 

The ghost propagator G is nothing but the inverse of the Faddeev-Popov operator M. 
We start our discussion in the continuum, where - in the case of covariant gauges - the 
Faddeev-Popov (FP) operator reads 

M = -V-D(A). (20) 

D(A) is the covariant derivative, emerging from the response of the gauge field to an 
infinitesimal gauge transformation. G is diagonal in the (adjoint) color space (a, b, . . . = 
1,...,8) 

G a \p) = 5 ab G(p 2 ) . (21) 
As an average over field configurations {^4}, it is translationally invariant 

G ab (x -y) = G(x - y) 5 ab = (M- 1 ^) ■ (22) 

In Landau gauge, because V • A = d^A^x) = 0, we can reverse the order of V and D, 

M — —V • D{A) — —D(A) ■ V . (23) 

To obtain the form of the FP operator on the lattice one closely mimics the steps one 
goes through in the continuum. The lattice form of that operator is nevertheless much more 
involved, since the response of the Lie-algebra A^ field to a gauge transformation does not 
have a closed form. In perturbation theory there is no logical obstruction, anyway: at any 
given order the lattice counterpart of the covariant derivative can be pinned down [2S1 ED] 
and all in all the FP operator reads 

M^=[d L ^V,(U)] ab 5 x , y . (24) 

Here we denote by T>(U) the lattice counterpart of the covariant derivative, indicating a 
dependence on U, even if this dependence is through a perturbative expansion. This results 
in a dependence on the adjoint representation of the Lie-algebra field 



1 + -$Jx) - —& u (x) - —$i(x) + — *— ®l(x) + 
2 M j 12 ^ v ; 720 mV ; 30240 M ; 



d* + i%(x). (25) 
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The $/j,(x) is the adjoint representation of the lattice gauge field <f> (<f> = — L4, A antiher- 
mitean) 

[%^)] bc =[tV^-if abc r,, (26) 
the fundamental representation being given by (i,j = 1, 2, 3): 

The ghost propagator in momentum space can in principle be obtained by Fourier 
transformation in each lattice configuration {U}, taking the translational invariance into 
account: 

G ab ( P (k)) = i / ^[M^g exp(ip(A;) ■ (x - y))\ = ([M" 1 ]" 6 > ( 28 ) 

\ x,y I {u} 

with the scalar product p(k) ■ x = ^ 2irk fJi X[ 1 / L^. Given the color-diagonal form, one is 
interested in the scalar function defined by the adjoint trace 

G{p{k)) = l -G aa {p{k)) = ~ (Tr adi M-\p(k))) {u} . (29) 

Due to its eight trivial zero eigenvalues (corresponding to global gauge transformations) 
the inverse M~ l needs to be defined with care. Since one is interested in non-zero momenta 
p(k), this means that one is automatically restricted to a subspace orthogonal to the space 
spanned by the (space-time constant) zero modes of the Faddeev-Popov operator. 

For an individual configuration, the inverse of the FP operator is not translation- 
ally invariant and hence ff28l is not so useful. The way out is quite standard and goes 
through the inversion of M on a delta source in momentum space, i.e. a plane wave source 
ijj(a, k) b x = 5 ab exp (i p(k) ■ x) in real space: 

[M-\k)) ah = ^rl>(a,k)?[M-%rl>{b,k) d y . (30) 

cd,xy 

As in standard non-perturbative applications, the issue is inverting the relevant operator 
on a source, but things are actually a lot easier in NSPT. The matrix M can simply be 
considered as expanded using the expansion of the gauge fields A">: 

M = + ^/T' /2 M« . (31) 
z>o 

Then the inverse of M will also obtained as a power series 

M" 1 = [M- l f ] + J2l 3 ~ l/2 l M ~ 1 f ) ■ (32) 

l>0 

One easily finds [20] a recursive way to obtain the higher orders of M~ l : 

[M _1 ] (0) = [M (0) r\ (33) 

i-i 

[M' 1 ]® = -[M- 1 ]^ J2 M(l ~ j) [M- 1 ]^. 

j=0 



It is worth pointing out that the computation proceeds back and forth from position to 
momentum space, since at any step the operators at hand are diagonal (or almost diagonal) 
in one of the two spaces. 

This is the time to reintroduce perturbative (loop) indices: we denote the n loop order 
of the ghost propagator as 

G^\ P {k)) = {[M-\ P {k))]f =2n) . (34) 

Thus a real loop order n of the ghost propagator corresponds to contributions with even 
powers (\//3) 2 , whereas contributions with odd powers {\ffi) 21+1 would indicate a non- 
loop contribution which does not exist by construction in standard perturbation theory. 
Contributions related to these non-integer n have to vanish numerically in NSPT for suf- 
ficiently large statistics. 

Similar to the gluon propagator in the forthcoming paper, we present the various orders 
of the ghost dressing function (or "form factor" ) in two forms: 

j(»)(p) = p 2 G {n \p(k)) , J(»>(p) = p 2 G {n \p{k)) . (35) 

The dressing function has to be calculated for each 4-momentum given by four integers - 
the momentum 4-tuple (hi, k 2 , k 3 , fc 4 ). 

Two comments are in order for ( |35l) . First of all, the J's are dimensionless quantities. 
Notice that in a numerical lattice simulation one does not have access to G(p(k)), but to 
the dimensionless quantity a~ 2 G(p(k)). Multiplying by (ap) 2 - or a 2 p 2 - one obtains the 
corresponding (still dimensionless) J(a,p) - or J(a,p) - we will be concerned with in the 
following. Second, the second notation in (135]) is chosen to remind that here the quantities 
p M are acting as the "physical" momenta. 

4 Relation to standard Lattice Perturbation Theory 

In the RI'-MOM scheme, the renormalized ghost dressing function J RI ' is defined as 

J [p,fi,a m/ ) = — r, (36) 

Z{a,(j,, am') 



with the condition 



J (p,H,a m >)\jp =l fl = 1. (37) 

Therefore, the ghost dressing function J(a,p, amO is j us t the ghost wave function renor- 
malization constant Z(a, //, am') a t — P 2 - 

The general form of an operator /function 0(a,fx,a) in terms of the renormalized cou- 
pling a(afj) = (g(an)) 2 / (16ti 2 ) @ can be written as 

i 

0(a, //, a) = 1 + a(a/i) 1 O it j (log(a/i))' ? . (38) 

i>0 j=0 



2 Note the additional factor l/(47r) compared to the "standard" definition. 
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Here j gives the possible powers of the logarithms in loop order i. 
Using (I36I) and (1371) . we get the expansions 

i 

Z(a } fi, a m ,) = l + Yl 4i' /Z z fj ( lo gM) j , (39) 

i>0 j=0 



J (a, p, a RV ) = 1 + ^ a nv z 3' i\ lo s(^ a ) 2 ) 

„'-n V Z ^ / 



(40) 



(41) 



j>0 j=0 

Since we measure the dressing function as an expansion in the lattice (inverse) coupling /3, 
we need the ghost wave function renormalization as an expansion in the lattice coupling 

i 

Z(a, fi, ocq) = 1 + ^ a o Y Zi 'i ( 1 °s( a / i )) J • ( 42 ) 

t>0 j=0 

From now we restrict ourselves to three-loop expressions in Landau gauge and in the 
quenched approximation. The logarithmic corrections are partly known from the anoma- 
lous dimension of the ghost field and the beta function and given as follows (compare e.g. 

3 N c , (43) 



2 



.hi' _ 35 2 m , _ 271 2 35 Rj/ 



z^=--N^ z^=~N^ + f^N e , (44) 



I..- 2765 „ T o Tiii 11933 „ T o 2765 i>i< . 
*S = - W N 3 C , zfl = N 3 C + — *JJ N 2 C , (45) 



RT , 157303 , „ 211 Ar . o 91 RT , Ar2 79 RT , 
*S = — iV c 3 + — C[3] iV c 3 + - z% Nl + — N c . 

The finite one-loop constant z^q has been calculated from the ghost self-energy in standard 
infinite volume LPT [32] with the value 

z«J = 13.8257. (46) 

To our knowledge, the constants z^o an d zfo' were not known so far. 

Taking into account the relation between the lattice coupling and the renormalized 
coupling in two-loop accuracy [331 EH [35] (the function d\ is scheme independent) 

am' = «ms + °( a Hs) = + d x ocl + d 2 a 3 , + • • • , (47) 



with 



di (a, = -2p log(a/x) + K x , (48) 
d 2 (a, fi) = (2/3 logM) 2 - 2(2/3 ^ + A) log(a/x) + K 2 , (49) 
f% = fN c , A = fiV c 2 , K 1 = -^ + 26.8384 iV c , 

K 2 = Kl + ^> - 452.047 + 197.252 iV c 2 , (50) 
10 



we get the coefficients of Z(a,//, «o) i n (H2|) up to three loops 



z l,l 


— z l,l 


) z l,0 — z l,0 , 










{01) 


z 2,2 


— z 2,2 


OR 7 Rl ' 7 - 7 Rl ' 4- ^ r RI ' 9/9 ^ RI ' 
— Zp Z x x , Z 2 ,l — Z 2X + Ai Z liX — ZfJ Z l Q , 








z 2,0 


- 7 RI ' _L K 7 RV 

— z 2 t- -f\l z 10 , 












z 3,3 


— z 3,3 














z 3.2 


- z RV 

— z 3,2 


+ 2K 1 z%-4[3 z% 


-4/3 #i 




- 2A 


z l,l ) 




Z 3,l 


- z m ' 

~ z 3,l 


+ 2Jjf! < + K 2 Z Z 


- 4A, *go 






z RV 


(53) 


z 3,0 


- 7 RI ' _J_ OK 7 RI ' 4- K 7 RV 

— Z 3,0 + /A 1 Z 2,0 + A 2 ^1,0 













Using these numbers and the inverse lattice coupling (3 we obtain the expression for the 
ghost dressing function to that order 

3 1 i 

J 3 - lo °P(a,p,/?) = 1 + £ p^Hw) , J (l) (a,p) = Yl J m H(pa) 2 Y , (54) 

i=l P j=0 

with the one-loop coefficients 

J X1 = -0.0854897, 

Ji,o = 0.0379954^ = 0.525314, (55) 

the two-loop coefficients 

J 2 , 2 = 0.0215195, 

J 21 = -0.313514 - 0.00324822 = -0.358423 , (56) 
J 2 = 0.106737 + 0.00144365 zjo = 1-47572 + 0.00144365 sgj , 

and the three-loop coefficients 

J 33 = -0.00660927, 

J 3i2 = 0.164130 + 0.000817642 zf^ = 0.175434 , 

J 31 = -1.38120 - 0.0210370 zgj - 0.000123418 

= -1.67205 - 0.000123418 zgj = -1.54589 - 0.0854897J 2 , , (57) 

J 3 = 0.375990 zgj + 0.00811105 zjj + 0.0000548523 zgj 

= 5.19833 + 0.00811105 zjj + 0.0000548523 z™j' . 

Thus we conclude: the coefficients of the leading logarithmic contributions for a given order 
can be exclusively taken from continuum computations whereas non-leading log coefficients 
are already influenced by finite lattice constants from corresponding lower loop orders. 

In general, these logarithmic terms are the relatively easy part of a diagrammatic 
perturbative computation, while getting the finite constants is the big effort. In NSPT 
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it is just the other way around. As it will be clear in Section 15. 2\ telling which is which 
in the measured ghost dressing function at each loop requires a fitting procedure: fitting 
a logarithmic term would actually require a terrific precision. As a consequence, we have 
just stated from the very beginning the (modest) attitude of not even trying to extract 
anomalous dimensions from our NSPT computations. On the contrary, we have to assume 
that the coefficients in front of the logarithms are known. Assuming that a continuum 
calculation has been performed up to a given loop order defines the maximum accuracy 
and the maximal loop number which we can reach in a lattice NSPT computation. 

We have to check that NSPT in the limits pa — > (or pa — > 0) and V — > oo reproduces 
the constant Ji$ from the one-loop result. After that we may proceed to determine from 
the two-loop measurements the unknown and, therefore, . Knowing that number, 
we are in the position to fix the coefficients of all non-leading logarithmic contributions 
in the three-loop expression and to estimate the constant J 3>0 or zf^ from the three-loop 
measurements. 

5 The practice of NSPT simulations 

5.1 Statistics and extrapolation to zero Langevin step 

Choosing a maximal loop expansion order n max = Z max /2 and solving the coupled system 
of equations for the gauge links Ux% and the gauge potentials A^ + ^i 2 simultaneously, one 
generates a sequence of these fields for all perturbative orders 1, . . . , l max . Such sequences 
- representing Langevin time trajectories — are generated for different finite Langevin 
time steps e and different lattices sizes, i.e. different N. 

From the perturbative fields at finite e we construct the observable of interest, in our 
case the ghost dressing function for each chosen set of lattice momenta and average it 
over the used configurations within a Langevin time trajectory. It is expected that the 
autocorrelation time for the observable is extending over subsequent configurations and 
increases with decreasing e. As a reasonable compromise between computer time and 
autocorrelation we have measured the ghost propagator using a sequence of configurations 
separated by 50 Langevin steps. Before measuring, Landau gauge has to be obtained on 
each configuration by satisfying (JTB1) to high precision. To further increase the statistics, 
we have generated several different trajectories for a given e and N. 

Since the ghost propagator has to be calculated for each individual momentum 4- 
tuple and the number of 4-tuples rapidly increases with the lattice size we have restricted 
ourselves to a certain set of momentum 4-tuples for our periodic symmetric lattice. We 
have chosen a cylindrical cut similar to Monte Carlo measurements such that the individual 
momentum components ap M (fc M ) do not differ too strongly from each other. The direction 
of that 4- vector ap is chosen nearly along a lattice 4-volume diagonal. In a 4-tuple we have 
chosen the ordered integers fc M modulo lattice size not to differ by more than one unit and 
have assumed a maximal difference |/c M — k u \ < 2. 

During the Langevin evolution we treated all 4-tuples leading to non-identical measure- 
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ments as independent ones. This led in our particular realization to 27 * (N/2) 4-tuples [f| 
for which we have measured the propagator. After obtaining the Langevin time averages 
we further averaged the propagator over 4-tuples which were equivalent due to lattice sym- 
metries and ended up with the following set of inequivalent momentum 4-tuples for our 
lattices chosen with even site numbers used in the further analysis: 

(yk . k j k . k^j . [k i kj k . k -0? ^ ^ 

N 

(k,k — l,k — l,k — 1), fc = l,...,-; (58) 

(k, k,k — l,k — 2), (Jfe, k - 1, k - 1, k - 2), 

N 

(k,k-l,k-2,k-2), k = 2,...,-. 

Choosing as example N = 8, out of originally chosen 108 4-tuples we get 25 inequivalent 
4-tuples for which the propagator has been measured. A final average has been performed 
on the results from different trajectories. 

In Tables [T] and [2] we present our collected statistics for different Langevin steps and 





N = 8 


N = 10 


N =12 


N = 14 


N =16 


N = 20 


e 


1 2 


1 2 


1 2 


1 2 


1 2 


1 2 


loop 


loop 


loop 


loop 


loop 


loop 


0.010 


2500 2000 


1800 1000 


1200 864 


1250 822 


800 384 


57 57 


0.020 


1750 1500 


1250 1000 


900 850 


750 650 


600 241 


43 43 


0.030 


1100 1000 


900 750 


750 778 


650 573 


373 262 


20 20 


0.050 


1100 1000 


900 900 


750 700 


650 650 


284 259 


39 39 


0.070 


1100 1000 


850 750 


700 770 


650 591 


259 192 


40 40 



Table 1: Number of one- loop and two- loop ghost propagator measurements for individual 
momentum 4-tuples (158]) at different lattices sizes done with the Leipzig NSPT code at 
Leipzig PC clusters. 

lattice volumes from the Leipzig and Parma code, respectively. Note that the Leipzig data 
contain only measurements up to two- loops (7 max = 4). Combining the measurements from 
all sources we have rather accurate data of the perturbative propagator at each momentum 
4-tuple. 

The behavior of the different orders of the dressing function J is shown in Figure [1] 
for loop numbers n = 1,2,3 and the non-loop contribution n = 3/2 vs. a 2 p 2 at different 
volumes and fixed Langevin step size e = 0.01: We note that the dressing function for 
the power P' 3 ^ 2 , which is not a loop contribution @, vanishes to a good accuracy for 
each individual lattice momentum 4-tuple after averaging over the configurations. This 
vanishing has to happen already at non-zero e. In general, the behavior in the infrared 

3 The described choice corresponds to the Leipzig code, in the Parma code not all of the highest momenta 
have been measured. 

4 By construction of the FP operator the contribution /3 -1 / 2 is zero in NSPT. 



13 



e 


AT 1 O 

iv — 1Z 


AT 1/1 

iv — 14 


AT if; 

iv — ID 


at on 

iv — ZU 


U.U1U, all loops 


yyu 


o/o 


1 11 A 

lzl4 




D D1 ^ nil lnnnc: 


n 

u 


n 

u 




R1 zL 


0.020, all loops 


753 


620 


968 


565 


0.030, all loops 


508 


438 








0.040, all loops 








539 


369 


0.050, all loops 


465 











0.070, all loops 


442 


301 









Table 2: Number of one-, two- and three-loop ghost propagator measurements for indi- 
vidual momentum 4-tuples fl58|) at different lattices sizes performed by the Parma code at 
Parma and Trento. 
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Figure 1: Measured ghost dressing function J vs. a 2 p 2 for all inequivalent lattice momen- 
tum 4-tuples near diagonal for N = 8, 10, 12, 14, 16, 20 and e = 0.01. Left: The one-loop 
two-loop (/3~ 2 ) and three-loop (P~ 3 ) contributions, right: the vanishing (oc (3~ 3 / 2 ) 
contribution. 

(small lattice momenta) is a bit more noisy. For a finite number of measurements the 
eventually still non-zero non-loop contribution n + | for some given 4-tuple is by several 
orders of magnitudes smaller compared to the "neighboring" n-th and (n + l)-th loop 
contributions. 

Now we are in the position to perform the extrapolation to zero Langevin step for 
the ghost dressing function at each individual momentum 4-tuple (ki, k 2 , k 3 , fc 4 ) from the 
available finite e measurements at fixed lattice size. Since the Langevin system is solved 
in the Euler scheme in which pieces of (e 2 ) are neglected, we fit the extrapolation of each 
order of the ghost dressing function to zero Langevin step by an ansatz of the form 
J(«) ( e ) = J(») ( p , L) + ae + be 2 . The fit result is a function of the lattice momentum 4-tuple 
and depends on the lattice size. For most of the cases we used 5 different e values, for the 
largest volume we had 7 Langevin steps for one- and two-loop and 4 such steps for the 
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three-loop extrapolation. 

The extrapolation to zero Langevin step is shown for momentum 4-tuples (2, 2, 1, 0) at 
TV = 12 and (1, 1, 1, 1) at N — 16 in Figures [2] and [31 respectively. In Figures [U we show 
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Figure 2: Extrapolation to e = of one-, two- and three- loop ghost dressing function 
for lattice size 12 4 and momentum 4-tuple (2,2, 1,0): the fitting function contains both a 
linear and a quadratic term in e. 




Figure 3: Same as in Figure [2] for lattice size 16 4 and momentum 4-tuple (1, 1, 1, 1). 
the extrapolation of the one- and two-loop contributions to zero Langevin step for N = 16. 

The remaining volume dependence of individual loops contributions vs. a 2 p 2 at zero 
Langevin step is shown in Figure for two-loop and three-loop dressing functions using 
some of the lattice volumes. 

One of the aims of this project is to provide referential perturbative results for finite 
lattices and for lattice momenta (k\, &2, k±) without any restriction. For this purpose, 
no further extrapolation is necessary. An exhaustive comparison with Monte Carlo results 
is not attempted in this paper. 

To reproduce LPT results in the infinite-volume limit at vanishing lattice spacing, 
however, different lattice sizes have to be studied in parallel as explained in the following. 

5.2 Fitting the data 

Our aim is to extract the finite constants J^o for loop order i in the perturbative expansion 
of the lattice ghost dressing function (l5"4"j) in infinite volume and in the continuum limit 
a — ► as accurate as possible. As pointed out in section |U at any loop, a computation 
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Figure 4: Extrapolation to zero Langevin step as a function of a p for a lattice of size 16 
left: one- loop dressing function, right: two-loop dressing function. 




Figure 5: Volume dependence at zero Langevin step as a function of a 2 p 2 for different 
lattice sizes; left: two-loop dressing function, right: three-loop dressing function. 

would yield a sum of a finite contribution plus logarithmic terms. To be definite, we recall 
that at one-loop order one gets 

J (1 W) = Ji,ilog(po) 2 + Ji, , (59) 

with the known coefficient J^i given in (155|) and similar expressions for higher loops. 
In general, at loop order % one gets terms up to power (log(pa) 2 )*, and all logarithmic 
contributions arise from anomalous dimensions and the beta function. Given what is 
known in the RI'-MOM scheme, in our present computation we can aim at three loops. 
In particular, reproducing the known one-loop lattice constant J\$ is an obvious check for 
our NSPT results. However, any NSPT computation is performed at both finite lattice 
spacing and finite lattice size. 

Let us first of all discuss the influence of the finite lattice spacing. Any lattice compu- 
tation does not have any explicit reference to the lattice spacing a; in a non-perturbative 
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computation this is always fixed a posteriori by matching to a physical scale. In the NSPT 
propagator computations one always handles dimensionless quantities which are fixed once 
a 4-tuple (ki, k<i, k 3 , k&) (with integer is selected. One has no direct access to a (dimen- 
sionful) momentum p^ (see ( TP7]) ). but only handles (notice the explicit reference to lattice 
size L = Na, with N a pure number) 

P,a = -jf • (60) 

While in the p^ definition the lattice size enters explicitly, one is left with the combination 
pa as variable, which for a ^ is the signature of finite lattice spacing also in infinite 
volume. Taking into account a 7^ 0, e.g. at one loop, the dressing function results in 

J (1) (pa) = J 1A \og(pa) 2 + J lfi (pa) , (61) 

instead of (1591) . The non-logarithmic constants Ji$(jpa) for loop number i are functions 
which encode the lattice artifacts. Since we expect these functions to be compliant to 
lattice symmetries, there is an obvious way to estimate the pa — > limit by means of a 
hypercubic-invariant Taylor series [36J. We collect measurements for (moderately) small 
values of pa and fit our results for the non- logarithmic part of a generic i-th loop to a 
formula which is fixed by the physical dimension and the scalar nature of the observable 
at hand 

(j)CL) 2 (j)(X) ^ 

Ji,o{pa) = Ji,o + Ci,i (pa) 2 + d,2 j—^ + cj.3 (pa) 4 + Q,4 {(pa) 2 ) + Q,5 j—^ H , (62) 

in which (pa) n stands for the hypercubic invariant ^^(Pm* 3 )™- This is a fit with 6 param- 
eters in which no (possible) subleading logarithms (proportional to 0(a 2 )) are taken into 
account 0. Apart from neglecting those logarithms, f|6^|) is a systematic expansion in pow- 
ers of pa. Going to next order, would require the inclusion of 5 extra parameters. Since J^o 
(i.e. the value at pa = 0) is the final number we are interested in, it is worth to stress that 
in principle the more we approach pa = 0, the better the fit is expected to work, at any 
fixed accuracy (in our case, O ((pa) 6 , . . .)). On the other side, the higher values of pa we 
take into account, the more terms in the Taylor series are expected to give a numerically 
significant contribution. 

Now we have to discuss how the finite lattice size has to be taken into account. In [36] 
it was pointed out that those effects can be large when an anomalous dimension is in place. 
Having at hand a variety of lattice extents, in this study we address a careful assessment 
of these effects (the main ideas entering the procedure can be found in [27] [22]). 

To take finite L effects into account, we consider as further generalization the ansatz 

J (1) (pa, pL) = J 1X log(pa) 2 + J lfi . L (pa, pL) , (63) 



J This is quite a strong hypothesis; the effectiveness of such a procedure can only validated a posteriori 
(e.g. by inspecting \ 2 values). In infinite volume these logs are calculable in principle. For the quark 
propagator and quark bilinears this has been done in [37j . 
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using the one-loop contribution as example. Dimensional arguments suggest a dependence 
on pL. Notice that this is not an irrelevant effect: a pL-effect would be there also in a 
continuum computation on a finite volume. In this notation Ji t o(pa) from (16T|) corresponds 
to the infinite lattice size limit J\fl(pa) = Ji^-^pa, oo) . Finite size effects can now be 
parametrized as 

J {1) (pa,pL) = J hl \og(pa) 2 + Ji, (pa) + (J lfi]L (pa,pL) - Ji,o(pa)) (64) 
= Ji 5 i log(pa) 2 + Jxjaipa) + 8J 1;0 (pa,pL) , 

which is, however, not yet useful for a fitting procedure. 

To make the ansatz usable for a fit purpose we once again perform the expansion (|62|) 
of J\fl{pa) in hypercubic invariants and take into account the (formal) continuum limit of 
5J lfi (pa,pL) 

5J lfi (pL) = 5J lfi (pa = 0,pL) . (65) 

Here pa corrections to the pL correction 8J\$(pL) are supposed to be corrections on cor- 
rections and will be neglected. 

For the non-logarithmic part of an i-th loop this procedure results in 

J ifi (pa,pL) = J ifi (pa) + 5J ifi (pL) + • • • , (66) 

where Ji t o(pa) is given by the expansion (1621 . Now fl66l) is indeed amenable to a fit. The 
key observation is the trivial equality ( 1601) p^L = p^aN = which can be treated as 
follows: Neglecting pa corrections to SJ it0 (pL) implies that measurements at the same fixed 
4-tuple (hi, &2, &3, k^) but at different lattice sizes L/a are affected by the same pL-effect 
(while pa values are different). 

This leads to a practical way of implementing the fit encoded in ( |66|) without assuming 
a functional form for the pL effects. Suppose we want to fit J\$(j)a,pL). First we select an 
interval [G»a)m in , (pa)^ ax ], say [0.25,6.25] for (pa) 2 . Then we select a collection of lattice 
sizes {iV = 8, 10, 12, 14, 16, 20} where measurements are available and identify all those 
4-tuples = (ki, k 2 , k 3 , k^) (resulting in values of (pa) 2 inside the chosen interval) which 
belong to each of the lattice sizes N = L/a. In the case at hand, these are 

{(1, 1, 1, 0), (1, 1, 1, 1), (2, 1, 0, 0), (2, 1, 1, 0), (2, 1, 1, 1), (2, 2, 1, 0), (2, 2, 1,1)}, (67) 

for which we choose the measured ghost dressing functions at one-loop after subtracting 
the logarithmic behavior J^i log(pa) 2 assumed to be known. For our 7 4-tuples and 6 
lattice sizes, we have just selected 42 measurements or data points. 

This choice so far still does not fix the overall normalization of the data. For that pur- 
pose we notice that finite- volume effects decrease both with increasing momentum squared 
and increasing lattice size. Therefore, we choose as reference fitting point - for a good 
approximation to pL = oo - a (log subtracted) data point from the largest available lat- 
tice at (pa) 2 approximately equal to the largest momentum squared of the chosen data 
set. In our example we add the measurement taken on the iV max = 20 lattice for the 4- 
tuple fc^max = (4,4,4,4). This choice leads to (pa) 2 = 6.31655 compared to the maximal 
(pa) 2 = 6.1685 from tuple (2, 2, 1, 1) at N — 8. All together we now have 43 data points 



18 



and 13 parameters to fit: 6 parameters J^o and {cxj ; / = 1, . . . , 5} and 7 5Ji t o(pL(k)) for 
each of the 7 4-tuples ^ & M max- By assumption we have SJi y o(pL ma ^(k mSuX )) = for the 
normalization point. 

The application of the described criteria to select measurements is analogous for higher 
loop cases. Before fitting the data in a chosen momentum range, all logarithmic pieces 
(supposed to be known) have to be subtracted. 

Finally let us discuss the effectiveness of our fitting procedure. A list of requirements 
to be fulfilled should be recalled: 

• A number of measurements on different N = L/a lattice sizes should be available. 

• An interval [(pa)^^, (pa) max] nas to be singled out in which a hypercubic Taylor 
expansion (with a manageable number of terms) for the quantity Jifl{pa) has to be 
effective. 

• The number of 4-tuples which enter our procedure has to be such that the total num- 
ber of data points is sufficiently large with respect to the number of fit parameters. 
That number has to be not too large. 

• It is assumed that the data point associated to /c max is virtually free of finite size 
effects. 

• All the procedure results in a sufficiently stable fit. 

The fulfillment of all these requirements has to be assessed a posteriori. Figures [6] and 
[7] shows the method at work at one loop and two loops. Looking at the left Figures l6|71 we 
observe that the numerical data (extrapolated to e — 0) from the chosen set of all lattice 
sizes - with all logarithmic contributions subtracted - scatter significantly. Switching off 
the finite volume 0{pL) corrections (SJifi(pL) in the fit form fl66l) ) in the right Figures I6HTI 
the stars (blue in color online) line up in 'rows' according to the different hypercubic 
invariants at infinite lattice volume. The reference point (here the rightmost point) is of 
course unchanged compared to its position in the left Figures. Finally, after removing also 
the non- rotational hypercubic artefacts (leaving only c^i ^ and c i4 ^ in (I66I) ). we 
obtain a smooth (almost linear) curve formed by the full circles (red in color online) which 
directly points towards the fitting constant J i)0 corresponding to the zero lattice spacing 
limit. The quality of the fitting procedure can be estimated by comparing the data with 
the fit points, as it can been seen in the left Figures. 

6 Results 

6.1 Finite constants in LPT up to three loops 

In order to determine the dressing function for the ghost propagator J(a,p, f3) up to three- 
loop order on the lattice we have to compute the constants Jk,o in the way outlined in 
Section 15.21 They fix all other coefficients of the subleading logarithms as it can be seen 
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Figure 6: Fitting of the one- loop non-logarithmic coefficient. Left: Data points after 
logarithmic subtraction for various lattice sizes N compared to points using the non-linear 
fit (1661) with (|62|) . Right: Stars are data after correction for finite volume effect, i.e. they 
represent Ji t o(pa); full circles are fit points after correcting both finite volume and some 
hypercubic effects with the exception of those proportional to the coefficients c^i and 
in fl66]). 
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Figure 7: Same as in Figure [6] for the two-loop non-logarithmic coefficient. 

in formulae ( 1551) . ( 1561) and ( 1571) . These relations also show that the coefficients of a certain 
order in perturbation theory depend on all lower-order results. 

One essential point in the fit procedure described in Section 15.21 is the existence of a 
stable fitting window (pa)mi n — (P a ) 2 — (P a )max- Practically, it is defined by the value of 
the quadratic error \ 2 of a nonlinear regression fit. It is not possible to give a prescription 
to fix a priori an optimal value for \ 2 besides the general condition x 2 > 1. It has to be 
determined for each case (i.e. fit function Jw(pa)) by inspection. In Figure [H] (left) we 
show the error x 2 as functions of (pa)^ in and (pa) 2 nax for the one-loop case. As discussed in 
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the example of Section 15721 the plot contains data for lattice sizes N = 8, 10, 12, 14, 16, 20. 
One recognizes several plateaus for the fit error x 2 - They prevent the determination of one 
optimal x 2 together with its attached fit parameters. But rather they suggest to average 
over the smallest possible x 2 ■ m Figure M (right) we show the values of the constant J^o 
for the smallest x 2 ^ 1 which agree within errors. 

Averaging over these data we obtain Ji j0 = 0.52523(95). We notice that this value is 
stable with respect to the introduction of weights proportional to the inverse of x 2 ■ The 
same holds for two and three loops. 

Another requirement for an efficient fit is a suitable relation of number of data points 
to the number of fit parameters. In Table [3] this can be read off for all considered orders. 
The value J\$ almost perfectly coincides with the exactly known constant given in (|55|) 




Figure 9: Same as in Figure [8] for the two-loop results, 
and the estimated error is small. This is a clear indication that we have reached a sufficient 
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Table 3: Number of fit parameters {Nf p ) and data points (Nd p ) for the x 2 taken to compute 
the constants Jk,o- 

good accuracy of the accumulated data in the NSPT approach and the proposed fitting 
procedure for handling both 0(pa) and 0(pL) effects works well. For two- loop order the 
corresponding results are shown in Figure [91 Averaging over the eight best \ 2 values we 
get J 2j0 = 1.4872(57). 

Three-loop data are available for lattice sizes L = 12, 14, 16, 20 only. Therefore, the fits 
have to be performed with a smaller set of input numbers. From Table [3] one can deduce 
that the three-loop fit has the least accuracy. The results are shown in Figure [101 The 
average yields Js t o = 4.94(27) which is supported by the Figure. 
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Figure 10: Same as in Figure [8] for the three-loop results. 

A major source of indetermination in our three loop analysis comes from the extrapo- 
lation e — > 0. For three sizes (N = 14, 16, 20) we only have four values of e. The number 
given above is obtained from a quadratic fit e — > used for one- and two-loop cases as 
well. Computing J^^ from a linear £-fit and using otherwise the same strategy we obtain 
5.40(11) which differs about 10% from the quadratic extrapolation. On the contrary, the 
one- and two-loop constants are practically not affected by the two extrapolations. For the 
linear extrapolation we get Ji >0 = 0.52451(74) and J 2)0 = 1.4867(43). 

Collecting all results we can write (1541) in a numerical form 

J(a,p,P) = 1 + - (-0.0854897 log(pa) 2 + 0.52523(95)) + 

P 

1 9 

+— (0.0215195 (log(pa) 2 ) - 0.358423 log(pa) 2 + 1.4872(57)) + (68) 
+— (-0.0066027 (log(pa) 2 ) 3 + 0.175434 (log(pa) 2 ) 2 - 1.6731 log(pa) 2 + 4.94(27)) . 

In (16 8 p the coefficients Jjj (j < i) are given with their significant digits only. 

Formulas fl56l) and fl57|) encode the relation to the and zfl coefficients. Due to 
numerical accidents, these are affected by relevant percental errors. 

6.2 A first comparison to MC data 

Using the A = —log U definition ([6]) as in NSPT in the gluon propagator and the Faddeev- 
Popov operator (1241) for the calculation of the ghost propagator, the Berlin Humboldt 
University group has obtained Monte Carlo results for the two propagators in Landau 
gauge and the resulting running gauge coupling [3B]. Maximizing the "linear" gauge fix- 
ing functional 8ltrU Xtfl as a preconditioner, the final gauge fixing was organized to 
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minimize 



Ew (*) E<£a 



(69) 



with an iterative scheme following (fTTj) . actually in a parallelized multigrid realization of 
Fourier acceleration (|T6|) . 

Since it is assumed that non-perturbative contributions dominate mainly the interme- 
diate and infrared momentum range, it is of interest here to compare directly the perturba- 
tive ghost dressing function obtained in NSPT with its Monte Carlo counterpart for each 
common momentum 4-tuple. 

We calculate the perturbative dressing function at a given lattice volume summed up 
to loop order n max for a given lattice coupling (3 as follows: 



(70) 



In Figure [TT] we compare the perturbative ghost dressing function at lattice size N — 16 
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Figure 11: Three- loop NSPT of the ghost dressing function in comparison to Monte Carlo 
at two (3 values. 

with Monte Carlo data at two different (3 values. We observe that at least three-loop 
accuracy is necessary to guarantee that the perturbative ghost propagator at larger f3 
((3 > 9.0) approximately describes the full two-point function in the large momentum 
squared region (p 2 > 1 GeV 2 ). 



7 Summary 

We presented two- and three-loop results for the ghost propagator of Lattice SU(3) in 
Landau gauge. A careful analysis of finite volume and finite lattice size effects has been 
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performed, whose methodology has been discussed quite in detail. The present work will 
be supplemented by a similar analysis for the gluon propagator (in the same gauge). This 
will open the way to a more careful comparison to Monte Carlo data for these quantities, 
which are supposed to encode informations on the confinement mechanism of non-Abelian 
gauge theories not only in the asymptotic infrared region, but also in the intermediate 
momentum range in the form of power corrections (condensates) [T3] and contributions 
from non-perturbative excitations (vortices etc.) [H],[T5] related to particular length scales. 
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